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Abstract: We investigate the four dimensional gravitational theories which admit 
homogeneous but anisotropic black brane solutions in asymptotically AdS space- 
time. The gravitational theories we consider are 1) Einstein-Maxwell-dilaton theory, 
and 2) Einstein-Maxwell-dilaton-axion theory with SL(2, R) symmetry. We obtain 
the solutions both analytically and numerically. Analytical solutions are obtained by 
perturbation from the isotropic solutions. Our solutions approach singular behavior 
at the horizon in the extremal limit but in non-extremal case, they are smooth 
everywhere. We also discuss how the third law of thermodynamics holds in our 
set-up. 
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1. Introduction 

Understanding the strongly coupled limit of quantum field theory is a long standing 
problem in theoretical physics, which exists in broad range of physics, from nuclear 
physics to condensed matter theory. The AdS/CFT correspondence, or more broadly, 
the gauge/gravity correspondence [1, 2, 3] is an extremely useful tool for this purpose 
since it gives a new perspective for strong coupled field theory from totally different 
theory viewpoint, i.e., gravitational theory on the asymptotically anti-de Sitter (AdS) 
viewpoint. Given this, it is very interesting to search for the generic gravitational 
solutions which are asymptotically AdS at UV but show non-trivial behavior at the 
interior, especially in IR. Each of these generic geometry corresponds to interesting 
phase in the field theory side. 

Famous such examples are black brane solutions in asymptotically AdS space- 
time. These solutions are very useful gravitational backgrounds for studying the 
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strongly coupled dual field theory, which corresponds to deconfined phase in QCD 
[4]. These black brane solutions, especially their charged ones, are also useful to study 
superfiuid phase [5,6,7] and "fractionallized" Fermi-liquid phase in condensed matter 
system [8], and they have been studied in great detail for the application to both 
QCD and condensed matter physics. For recent review of these, see for examples, 
[9, 10, 11, 12, 13, 14, 15, 16]. 

Generically, finding generic black brane solutions of Einstein equations are quite 
difficult task, therefore we usually put several ansatz for the metric to simplify equa- 
tions of motion. As a result, we end up with studying holographic model with highly 
symmetric restrictions on the field theory side. For example, one of the well studied 
system takes the metric form such as 

ds^ = -a^(r)dt^ + — — + l?ir)idx^ + dy'^) , (1.1) 

which corresponds to the homogeneous and isotropic system in the dual field theory, 
and one of such solutions is, of course, well-known Schwarzschild black brane solution 
or Reissner- Nordstrom solution in AdS. However in more generic setting like Einstein- 
Maxwell-dilaton gravity, we can obtain more generic values for a(r) and h{r). One 
such nice example is the Lifshitz geometry, which has a scaling symmetry but this 
scaling symmetry acts nontrivially. See, for example for the application of such 
generic Lifshitz-like geometries from Einstein-Maxwell-dilaton gravity to condensed 
matter system, [17] - [32]. 

However in the condensed matter systems in nature, there are systems which 
do not satisfy above homogeneity and isotropy. For example, due to the crystal 
structures of the background atomic lattices, some systems induce the anisotropic 
structure of the Fermi-surfaces for the electrons. Generically, Fermi surfaces are 
generically not spherically symmetric but rather quite anisotropic in momentum 
space. Given that such anisotropic systems are quite ubiquitous in condensed matter 
theory, obviously it is desirable to study more generic gravity solutions where we can 
relax above homogeneity and isotropy condition at IR, and study such systems for the 
holographic setting. In this paper we would like to study more gravitational solutions 
which are homogeneous but not isotropic at the fixed radial slice. Putting asides the 
holographic applications, studying the generic solutions in the gravity system, which 
are homogeneous but not isotropic, is, by itself, very interesting problems. 

Recently, by applying the Bianchi classification well-studied in cosmology, it is 
shown in [33] that we can obtain very generic new classes of solutions which are 
homogeneous but anisotropic in IR. In this paper, we restrict our attention to the 
following four-dimensional metric ansatz; 

ds"" = (-g{z)de + ^ + e^(^)+^(^)dx2 + 2c{z)dxdy + e^(^)-^(^)rfy2^ ^ (^ 2) 
\ g[z) J 
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This is the generahzation of the ansatz (1.1)^, and admits static, homogeneous, but 
anisotropic metric where homogeneity at the fixed radial shce is achieved through two 
Kilhng vectors dx-, and dy\ The new function B{z) and c{z) are introduced to allow 
anisotropy. The point that all functions are dependent on only radial coordinate z, 
guarantees that this system is homogeneous at the fixed radial slice, therefore it is 
so at the dual field theory side.^ Note that we can make the metric diagonal in x, y 
place only at some fixed radius z = Zi. This means that by coordinate transformation 
for X and we can make B{zi) = c{zi) = only at some radius z = Zi. However, 
for generic metric and generic radius, we cannot make the metric diagonal in (x, y) 
place and this induces the anisotropy to the system. 

The organization of this paper is followings; In section 2, we start with the anal- 
ysis of what kinds of matter profiles allow anisotropy for the homogeneous metric 
ansatz. Quite obviously by adding more new degrees of freedom with more generic 
profile, we can always construct more complicated anisotropic solutions. As many of 
the bottom up holographic approach does not have a definite principle for what de- 
grees of freedom are needed, we consider the case where anisotropy is simply induced 
by the scalar fields. In section 3, we will construct both perturbatively and numeri- 
cally some of such anisotropic solutions in the Einstein-Maxwell-dilaton system. In 
section 4, we generalize our construction to the Eistein-Maxwell-axion-dilaton sys- 
tem which allows SL{2,R) symmetry. The good point is that once we obtain one 
solution, by using the SL{2,R) symmetry, the symmetry allows us to construct an- 
other solutions. Section 5 is devoted to the argument for horizon area under the 
presence of anisotropic matter. We show how the thermodynamical third law holds 
in our set-up, implying that the horizon area goes to vanish at low temperature. We 
end with discussion at section 6.^ 

2. Searching gravitational theories admitting anisotropy 

We start with the following question; from what kind of gravitational theories, can 
we have homogeneous but anisotropic black brane solutions on the metric form (1.2)? 
Obviously by adding more and more degrees of freedom, we can have more and more 

^Here we have rescaled the radial coordinate z as 1/z = r, such that we have exphcit overall 
factor l/z^, just for convenience. 

^We can obtain homogeneous space-time without two commutative Killing vectors dx, and dy. 
In more generic situations, we can consider Bianchi types of the geometries where homogeneity 
is achieved by two non-commutative Killing vectors. See [33] for the five-dimensional analysis, 
where three spatial Killing vectors do not commute, except for the type I case. In this paper, we 
consider the simplest case, where metric admits two commutative Killing vectors dx, and dy. This 
is analogous to the type I case studied in [33]. 

^For another approach for anisotropic black brane study, see [34, 35]. See also [36, 37, 38] 
for spontaneous homogeneity symmetry breaking in holographic approach. Recently holographic 
spontaneous isotropy breaking in M-theory setting is discussed in [39]. 
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complicated anisotropic solutions. Therefore we first ask what is the simplest grav- 
itational theory admitting such anisotropic solutions. For concreteness, we restrict 
out attention to the 3 + 1 dimensional static homogeneous black brane solutions in 
asymptotically AdS spacetime. 

2.1 Generalized Gaussian Null Coordinates 

Before we analyze each system in detail, we first comment on the metric ansatz we 
choose. For the static regular homogeneous black branes, we can make the metric in 
the form written by a generalized gaussian null coordinates 

ds'^ = — g^y dx'^dx" 

= ^ i^g{z)dv^ - 2dvdz + e^(")+^(")c/x2 + e^^-'^'^^^Uy^ + 2c{z)dxdy^ , (2.1) 

where B[z) or c{z) generates the anisotropy of the metric. This metric can be brought 
back to the form (1.2) by coordinate transformation 

dz 

dv = dt , 2.2 

where g{z) > outside the horizon. 

On the generalized gaussian null coordinate, the horizon is represented as 2; = 
constant null hypersurface where g{z) = 0. Since the determinant of the metric is 
non-zero at the horizon, the coordinate form has advantage over the horizon such 
that there is no coordinate singularity on the horizon. In addition, the vector field 
d/dz is the generator of a null geodesic curve and z is the afiine parameter on the 
conformal metric g^^. So, it is natural to assume that all the metric functions g{z), 
A{z), B{z), and c{z) are with respect to z. As shown in Appendix B, all the 
scalar curvature squares are finite if all the metric functions are and e^"^ — > 0. 
The latter means that the physical area spanned by x and y are finite. Hereafter, 
without loss of generality, we can set the location of the horizon to z = 1, then 

g{z = 1) = 0. (2.3) 

2.2 Pure Gravity 

We first consider the pure gravity system with action 

S= f d'x^g f + 4 V (2.4) 



and ask if this theory admits anisotropic solution of the type we mention in the 
introduction. Here L is the AdS curvature radius. 



-4 - 



With the coordinates (2.1), the Einstein equations of pure gravity system (2.4) 
can be written, after the several massage, as follows 

g2A^2^g2A _ ^2^gj^, ^ e^^z[g{-Ae^^ + lOc^ + zee') + ^(e^^ - 3e^)g']A' 

+e^^{g{6 + z^A'^) - 2zg' - 6) + e^[z^gc^ + 2zcc\zg' - Ag) - 40^(3 - 3^ + zg')] 

+e^^[2zcc'g - z^gc'^ - 2c'^{-9 + g{9 + z^A'^) - 3zg')] = 0, (2.5) 

z{c' - e'^)gB" + [e'\g{2 - zA') - zg') 

-c{zgc' + ci2gizA' - 1) - zg'))]B' = 0, (2.6) 

z\e^^ - c'')gc" - z[g{Ac^ + €^^(2 + zA)) - z{^^ + c^)(^Y 

+ [z^gd' + c\-Q + &g- 2zg') + 26^^(3 + {zA - l){?,g - zg'))]c = 0, (2.7) 

z'^gc'^ + 2zcc'{zg' - Ag) - 4c^(3 - 3^ + zg') 

+e^^[12 + g{-12 + z{8A' - zA'^ + zB'""^)) + 2^(2 - zA')g'] = 0, (2.8) 

where the last equation corresponds to the constraint equation. 

Let us consider initial data a.t z = Zi timelike hypersurface (0 < Zi < 1). By 
using the coordinate transformation in [x, y) plane, we can make the metric aX z = Zi 
such that 

B{z,) = c{zi) = Q, (2.9) 

and then, {x^y) plane metric becomes ds\ = e^^^^\dx'^ + dy"^) at 2; = Zi. Furthermore, 
by an orthogonal transformation at z = Zi, we can rotate coordinates {x,y) so that 
we can set c'{zi) = 0. So, initial data at 2; = 2;j timelike hypersurface is reduced to 

A{z,), A'izi), B'izi), gizi), B{z,) = c(z,) = c'{z,) = . (2.10) 

Since e^^^^^ - c'^{z) > and g{z) > for < ^ < 1, the Eqs (2.5) - (2.7) are 
the regular second order differential equations. So, given the initial data (2.10), the 
solutions of the Einstein equations are uniquely determined for < 2 < 1. Now 
suppose that c{z) = for < 2; < 1. Then, substituting c{z) = into Eq. (2.6) and 
integrating once by z, one obtains 

Q 2 -A{z) 

B'i^)-^^, (2.11) 

where Ci is a constant of integration. By assuming regularity of the metric at the 
horizon 2 = 1, we must set Ci to zero, as g{l) = 0. Therefore, with (2.9), we will 
obtain a solution with B{z) = c{z) = 0, which implies isotropy. However due to the 
uniqueness of the solution, this is the unique solution with a regular event horizon of 
Eqs. (2.5) - (2.8) satisfying the initial data (2.10). This means that all the solutions 
with a regular event horizon must be homogeneous and isotropic in the pure gravity 
system ^. 

"^We can say that the B{z) = c{z) = solution is also unique just inside the event horizon. We 
give the argument in Appendix A. 
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To see there is a homogeneous and isotropic solution satisfying initial condition 
(2.10) with B{z) = c{z) = in more detail, let us substitute B{z) = c{z) = 0, 
and then Einstein Eqs. (2.5) and (2.8) are reduced to the following two coupled 
differential equations 

z^gA" + z{zg' - Ag)A' - 6 + ^(6 + z'^A'^) - 2zg' = 0, (2.12) 
12 + g{z{8A' - zA'^) - 12) + 2z{2 - zA')g' = 0, (2.13) 

while Eqs. (2.6) and (2.7) are trivially satisfied. Under the coordinate transformation 
in {x,y) plane at 2; = 1, we can set A{1) = 0. From Eqs. (2.12) x2+(2.13), we obtain 

2A" + A'^ = 0. (2.14) 

The solution satisfying A(l) = is 



A 



A'{1) [ A{1 



2 



2+1-^ . (2.15) 



Putting back this into (2.13) and by choosing the integration constant such that 
5^(1) = is satisfied, we obtain 

_ Ajlfiz - If + 2A'{1) jz' -3z + 2)+Az'-A 

(A'(l)-2)(A'(l)(^-l) + 2) ■ ^ ^ 

Now let us consider the asymptotic AdS boundary condition. At 2; — 0, we have 



2A'(l)z A'ilfz^ 8z^ ^, , 

Therefore, putting the asymptotically AdS boundary condition at z = requires 
A\l) = 0. Therefore from (2.15), we have A{z) = 0. Substituting A'{1) = into 
Eq. (2.16), g{z) becomes 

g{z) = l-z\ (2.18) 

Introducing time coordinate t as (2.2), we obtain the familiar form of Schwarzschild- 
AdS metric, 

ds' = ^(^-{1- z')dt' + + dx' + dy']. (2.19) 

Thus, in this pure gravity set-up, due to the fact that the bulk degrees of freedom 
is tiny, it is not possible to obtain the regular solution whose horizon is regular and 
homogeneous but show anisotropy as ansatz (2.1). Therefore we consider adding 
matter degrees of freedom and see if we can relax this constraint. 
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2.3 Einstein Equations with General Matter Fields 

Next we consider the system consisting of gravity and generic matter with the action: 

S = J d'^xV^ (^R + + CU9,dg)y (2.20) 

where Cm are the Lagrangian density of the generic matter. The evolution equation 
for c{z) in the metric (1.2) is derived from the Einstein equations as: 

^2(g2^ - c^)gc" - z[g{Ac^ + e^^{2 + zJ^)) - z{^^ + c^)g'Y 
+ [z^gc'^ + c^{-6 + 6g- 2zg') + 26^^(3 + {zA' - l){3g - zg'))]c 
+L\e^^ - c2)(3T^ - T\ + T\ - Tyy)c = -2L\e^^ - c^)T\ , (2.21) 

where T^^, = —5Cm{g,dg)/5g^'' is the energy-momentum tensor. This equation tells 
us that c{z) is always generated by the non-zero component T^^. Suppose that 
T^x = 0. Then, we can obtain solutions with c{z) = 0. For simplicity, hereafter, we 
shall restrict our search for solutions which satisfy c{z) = in this paper. ^ 
Then, substitution of c{z) = into the Einstein equations restricts T^^ as 

T^^ = for any fi u except T\ . (2.22) 

The other non-trivial Einstein equations are reduced to, after several massage, 

2z V + z^9{A'^ + B'^) = 2L\T\ - T\) , (2.23) 
z^gB" + z{zg' + {zA' - 2)g}B' = L^{T\ - T\) , (2.24) 
2z{zA' - 2)g' + 12(^ - 1) + zg{A!{zA! - 8) - zB'^} = AL^T\ , (2.25) 

where the last equation corresponds to the constraint equation. The second Eq. (2.24) 
implies that non-zero B{z) can be generated by the anisotropic energy momentum 
tensor in which 

T\^T\. (2.26) 

In this paper, we seek the gravitational system with matter satisfying (2.26). 

Let us consider, for example, the system consisting of a real scalar field and 
gravity with action: 

S = j d'xV^ {^^J2- - ^('^)) ' (2-27) 

^Actually, as shown in section 2.2, wc can show that c(z) = is the unique solution satisfying 
the initial data c(zi) = c'(z,) = in Eq. (2.10) when T^.^ = 0. 
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where V{(f)) is the potential of the scalar field (p. The field Eqs. of the action (2.27) 
are 

1 3 

T^u = V^0V,0 - ^g^, ((V0)2 + r(0)) , (2.28) 
□0 = ^V'{<P). (2.29) 

Suppose depends only on z, then we have T^x = T^y and in quite analogous 
to the pure gravity case, we can conclude that B{z) must be zero. This implies that 
there is no anisotropic solution due to the uniqueness of the solution. 

In the next section, we consider the matter fields such that (2.26) is satisfied and 
can have homogeneous but anisotropic solutions. 



3. Anisotropic black branes in Einstein-Maxwell-dilaton the- 
ory 

In this section, we analytically and numerically obtain anisotropic black brane so- 
lutions in Einstein-Maxwell theory coupled to a real massless scalar (dilaton) field 
0. We introduce gauge potential in the bulk so that it plays the role of chemical 
potential in the dual field theory. The action we consider is; 



S 



d'^x^J—g 



R + 



6 



(V0) 



2 p pf^i^ 



with F^i, = df^Ay — dyA^. Field equations are 

R. 



□0 = 0, 
^/^^V^F^'^ = d^iy^F^") = . 

We shall take an ansatz for the metric (2.1) 

L2 



(3.1) 



(3.2) 

(3.3) 
(3.4) 



ds 



-g{z)dv^ - 2dvdz + e^^'^+^^'Ux^ + e^^'^-^^'^^dy^ + 2c{z)dxdy , (3.5) 



and the gauge field as 



A^dx^ = Ay{z)dv 



(3.6) 



Just for simplicity, we furthermore restrict our attention only to the solutions which 
satisfy 



ciz 



0. 



(3.7) 
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As shown in section 2, we find that the system cannot have anisotropic regular 
black brane solutions when (j) has only radial coordinate dependence as = (piz). To 
see this, note that the energy-momentum tensor for the matter field from the action 
(3.1) is 

n = ^^(/.V.*/. - ^5^(V</))2 + If\F,' - . (3.8) 

Given that gauge potential is a function of radial coordinate, we can have nonzero 
T^x — T^y if scalar field has coordinate x- or ?/-dependence. If we choose for the 
scalar field as [34, 35], 

(p = ax (3.9) 

where a is a constant, then 

T\-T% = aV^7^0. (3.10) 

Of course, the profile (3.9) implies T^^ = 0. Therefore it is consistent with the 
assumption (3.7). So, we consider a ^ for anisotropic solution ^. Note that we 
have introduced the manifest x-dependence for the scalar field 0, but not for the 
metric and gauge potential. Therefore, metric is still homogeneous, even though the 
scalar field induces inhomogeneity. 

Under the ansatz, Eq. (3.3) is automatically satisfied. The solution of Eq. (3.4) 
becomes 

= A[, = Lcii^ig"'')^)-' = Lcie-^, (3.11) 

where ci is a constant corresponding to charge density of the gauge field. Substituting 

the ansatz (3.6) and Eq. (3.11) into the Einstein Eqs. (3.2), we obtain the following 
three coupled differential equations, 

2A" + A'^ + B'"^ = 0, (3.12) 

zgB" + {zg' + g{zA' - 2)}B' + a^ze"^^^ = 0, (3.13) 
2z{zA' - 2)g' + {12 + zA'{zA' - 8) - z^B'^}g 

+ clz'e-^^ + 2{a^z'e-^~^ - 6) = 0. (3.14) 

Since we seek the black brane solutions with smooth horizon, we require all the 
metric functions are C^. Before we investigate these equations of motion to find 

^More explicitly, for anisotropy we need a ^ can be seen as follows; by rescaling coordinate y 
as y ^ ay, B{1) can be set to zero. By Eq. (3.13), we also obtain B'{1) = 0. Then, we have initial 
condition where only one parameter is given by ^'(1), while g(l) = and g'(l) is determined by 
(3.14). This initial condition uniquely determines the solution, and it is straightforward to obtain 
solution with the ansatz B{z) = 0, which is Reissner-Nordstrom solution in AdS. The regular 
solution obtained with B{z) — obviously satisfies the boundary conditions B{1) = B'{1) = 0, 
therefore they are the unique regular solutions of Eqs. (3.12) - (3.14). 
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explicit solutions, we point out that there is a difference between non-extremal case 
and extremal case of the system. For a 7^ 0, the above equations of motion at the 
horizon {z = 1), which we call "horizon condition" becomes 

2A"{1) + A'iiy + B'{lf = (3.15) 
g'{l)B'{l) + a2e"^(i)-^(i) = (3.16) 
2{A'{1) - 2)g'{l) + c?e-2^« + 2(a2e-^W-^W - 6) = (3.17) 

with the assumption that e"^'-^-', e^*-^^ are non-singular. This "horizon condition" 
tells the crucial difference between extremal limit (zero temperature limit) and non- 
extremal case (non-zero temperature). If black branes are extremal, then around the 
horizon z = 1, g{z) is expected to have double zero, as 5^ ~ (1—2;)^. Then, g' ~ (l—z), 
and from (3.16), we see that B'{1) diverges as 1/ (1 — 2), so -B(l) diverges as log(l — 2;). 
Note that these divergence contradicts with our non-singular assumption for -B(l), 
and this strongly supports the anisotropy diverges at the extremal limit. On the 
other hand, for non-extremal case, g ~ {1 — z), and g' ~ const. Therefore -B'(l) 
approaches some constant, which is consistent with the non-singular assumption for 
B{1). This crucial difference between non-extremal case and extremal case implies 
that extremal limit might be singular. 

3.1 Analytic solutions by perturbations 

When a is very small, we can obtain the analytic solutions of Eqs. (3.12) - (3.14) by 
perturbative expansion in a. The unperturbed isotropic black brane solution is the 
Reissner-Nordstrom solution in AdS, given by 

ds^ = ^ \ —gQ{z)dv^ — 2dvdz + dx^ + dy^ 

go{z) = l-{l + 'fjz' + 'lz\ (3.18) 
For later convenience, let us define the non-extremal parameter ^ as 

ci = A{i+e+e)- (3.19) 

Then, the function g^ is rewritten by 

g,{z) = 1 - (1 + e + + ey + {i+e+ e)z' 

= {l-z){l- iz){l + (1 + e)2 + (1 + e + e)z^}. (3.20) 

where the event horizon and the inner horizon are located as 2; = 1 and z = 1/^, 
respectively. Therefore, < < 1 and the extremal limit corresponds to = 1 (c^ = 
12). 
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Let us expand the anisotropic solutions of Eqs. (3.12) - (3.14) around the Reissner- 
Nordstrom solution in AdS as a series expansion in a: 

A{z) = Ao{z) + a'^Ai{z) + 0^^2(2) + ■ ■ ■ , 
B{z) = Bo{z) + a^Bi{z) + a^B2{z) + ■ ■ • , 
g{z) = go{z) + a^g^iz) + a^g^^z) + ■ ■ ■ , (3.21) 

where Aq{z) = Bq{z) = 0. Substituting Eqs. (3.21) into Eqs. (3.12) - (3.14), we 
obtain equations of motion for Ai, Bi, gi, as 

A'l = 0, (3.22) 

zgoB'l + {zg'^ - 2go)B[ + z = 0, (3.23) 

g[--^gi-'- = 0. (3.24) 

The regular solutions of Eqs. (3.22), (3.23), and (3.24) are easily obtained as 

z'^ 

Ai{z) = ca,,z + ca,, , gi{z) = - — {1 -Cg^z) , (3.25) 

0/ ^ z{cb,z + 1) 

' {l-z){l-^z)iz{l+^ + z + ^z + ez) + l)' ^' ' 

where ca^^^ ^a^j, c^i, Cg^ are constants to be determined from the boundary condition. 
By perturbation expansion in a^, (3.16) yields 

g[{l)B',{l) + g',{l)B[{l) + e-^o«-^«(i) = . (3.27) 

Using Ao(l) = 5o(l) = 5^(1) = 0, and ^^(1) = -3 + ^ + + e^ we have 

^i(l) = 3_^_\2_^3 > (3-28) 

for ,^ < 1 case. This condition determines that cb^ in (3.26) must be 

cb, = -1, (3.29) 

otherwise, Bi[z) diverge as log(l — z) at the horizon z = 1. 

In the extremal case ^ = 1, we face the breakdown of perturbation as we pointed 
out before. In that case, 

gextremal^^^ = (1 - 2)^(1 + 2z + 3z^) , (3.30) 

therefore, g^il) goes to zero as 1— Then B[{z) diverges as 1/(1— 2;), therefore, -Bi(l) 
diverges logarithmically as log(l — z). However this contradicts with the assumption 
that all functions are smooth at the horizon. Since perturbation breaks down at the 
extremal limit, we will consider the non-extremal case (^ < 1) from now on. 
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Since z = 1 is horizon, we need gi{l) = 0. This imphes that in (3.25), 

c,, = 1 , g,(z) = -'-^{l-z), (3.31) 

Then, g[{l) = 1/2. Perturbation expansion in of (3.17) with Aq{z) = Bq{z) = 
gives 

2A[{l)g',{l) - Ag[{l) - 2clA^{l) + 2 = 0. (3.32) 
With g[{l) = 1/2, (3.32) becomes, 

e)cA,^ = 4(e +e+ e)cA,, . (3.33) 

Finally we impose boundary condition at 2; = that the solution is asymptotic to 
AdS. We can require the boundary condition such that 

Ai(0) = A;(0) = 0, (3.34) 

which is consistent with (3.33). This ends up determining all the constant and we 
finally obtain, 

A,{z) = 0, g^{z) = -^{l-z), 

2 ln(l - ^z) - ln{l + (1 + e)^ + (1 + e + e)^'} 

-DlU) 



2(1 + 2^ + 3^) 
1 + 3^ 



(3.35) 



(i + 2e + 3e)v/3 + 2e + 3e 

f arctan ( , ] - arctan ( ^ ^ ^ + + ^ + e)^' 

This satisfies asymptotic AdS condition -Bi(O) = B[{0) = as 

B^iz) = ^ + 0{z^) . (z^O) (3.36) 

Similarly, we can easily show that the second order solutions, A2, B2, and g2 
are also regular outside and on the horizon as follows. Substituting Eqs. (3.21) into 
Eqs. (3.12) - (3.14), we obtain equations of motion for A2, g2, and B2 as 

A'i = ~, (3.37) 

, 3 zAW. f „ zB\^\ zBi ^ , , 

92 - -92 = - + ^ J g^-^^S,, (3.38) 

+ f- - -V2 = - + f — - -Vi - -^1 = ^B. (3.39) 
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By imposing ^2(0) = ^2(0) = at the boundary condition for asymptotically 
AdS spacetime, we can formally obtain the solution A2 by Eq. (3.37) as 

(3.40) 

which is obviously regular outside and on the horizon, z = 1, a.s B[ is finite there. 
This solution means that the horizon area per unit coordinate interval, Ax = Ay = 1, 
e^(i) decreases due to the existence of anisotropy, B[ ^ 0. 

The solution of Eq. (3.38) is also obtained by the boundary condition (72(1) = 

as 

g.{z) = r ^dz'. (3.41) 
Ji ^ 

Since Sg is regular outside and on the horizon, and decays as Sg ~ z"^ near the 
boundary of asymptotically AdS spacetime, §2 is also regular there. Finally, B2 is 
also obtained by integrating Eq. (3.39) as 

where we imposed a regular boundary condition on the horizon. Using the fact that 
Sb is regular outside and on the horizon and it decays as Sb ~ near the AdS 
boundary, we find that B'2 is also regular there and -82(0) = 0. B2 is obtained by 
imposing the boundary condition i?2(0) = as 

Thus, we have checked that the perturbed solutions in Eqs. (3.21) expanded in a 
series of satisfy the boundary condition near the AdS boundary: 

A{z) = 0(^2), B{z) = 0(^2), g{z) = 1 + 0{z'), (3.44) 

and they are regular outside and on the horizon, up to the second order O(a^). 

We can in principle continue working on higher corrections and obtain anisotropic 
solutions by perturbation in this way. 

3.2 Numerical solutions 

We can also numerically obtain anisotropic black brane solutions by solving Eqs. (3. 12) 
- (3.14). By rescaling x and y properly, A{1) and -B(l) can be set to zero. Under 
the conditions, we can read off the regularity condition at the horizon from "horizon 
condition" Eqs. (3.15) - (3.17) as 

O'^ — fir? o*^ 
A'il) = 2+- - + ^, B' {!) = -, (3.45) 
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where k, is defined as k = —g'{l). Thus, the regular solutions of Eqs. (3.12) - (3.14) 
are uniquely determined by three parameters, a, ci, and n. 

As the boundary conditions for the asymptotically AdS spacetime, we require 

that 

A{0) = B{0) = A'{0) = B'{0) = . (3.46) 

Under the boundary conditions, one obtains the asymptotic behavior of g{z) by 
Eq. (3.14) as 

g{z) = 1 + 0{z^) . (3.47) 

The first two conditions of Eqs. (3.46) are automatically satisfied by rescaling x 
and y again for the solutions of Eqs. (3.12) - (3.14) satisfying the regularity condi- 
tions (3.45), 

^ ^ g-(A(0)+B(0))/2^ ^ y ^ g-(A(0)-B(0))/2^ ^ ^ ^ ^^_(A(0)+B(0))/2 _ ^3_4g) 

So, the only task is to obtain the solutions satisfying the latter two conditions, 
A'{0) = -B'(O) = 0. -B'(O) = is automatically satisfied because the asymptotic 
behavior of B is derived from Eq. (3.13) as 

B{z) = B{0) + 0{z^) . (3.49) 

Thus, the anisotropic black brane solutions in asymptotically AdS spacetime are 
obtained by searching the parameters a, Ci, and k satisfying A'{1) = 0. 

We numerically find the value ci satisfying A'{0) = for a fixed a and k by 
solving Eqs. (3.12) - (3.14) from the horizon {z = 1) to the infinity {z = 0) under 
the regularity condition (3.45). Thus, the numerical solutions rescaled by Eq. (3.48) 
always satisfy the boundary conditions (3.46). 

Figs. 1-6 show the anisotropic black brane solutions for the rescaled coordinates. 
Figs. 1-2 plot the horizon area e^^^^ and -B(l) as a function of k, (the minimum of k, 
is 2 X 10"^) for a = V2. As shown in Fig. 3, the curve B{1) = 0.862 fits the 

plot of Fig. 3 well. This suggests that -B(l) diverges as k, goes to zero. According 
to the divergence, anisotropy grows near the horizon as the solution approaches 
extremal (k — )■ 0). Figs. 4-6 show the metric as a function of z for each k and 
a = V2. 

By Eqs. (3.12) - (3.14), we can argue that the horizon area e^^^^ per unit length. 
Ax = Ay = 1 must go to zero as k — > 0. Eq. (3.12) is rewritten in the form of 
Raychaudhuri equation for the null geodesies: 

^' = -V- V, (3.50) 
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Figure 1: (color online) K-e^^^^ relation. The area of the horizon per unit Ax = Ay = 1 
are shown function k for a = \f2. 
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Figure 2: (color online) k-B{1) relation for a = \/2. B{1) rapidly grows as k decreases 
to zero. 
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Figure 3: (color online) k-B{1) relation near k = for a = \f2. The solid curve B(l) = 
0.862 fits the plot well. 



where Q is the expansion of the congruence defined as 6* = A'. Integrating this by z 



D 




0.2 0.4 0.6 0.8 l.( 



Figure 4: (color online) D = z^gxT/L"^ (blue), z^g.yy/L?' (green), and g (red) are shown 
as a function of z for a = \/2, ^ = 0.5. 
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Figure 5: (color online) D = z'^gxx/L'^ (blue), z'^gyy/L?' (green), and g (red) are shown 
as a function of z for a = \/2, n = 0.1. 
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Figure 6: (color online) D = z'^gxx/L'^ (blue), z'^gyy/L?' (green), and g (red) are shown 
as a function of z for a = \/2, k, = 0.001. 



twice, we obtain the inequality 

A{z) < - 




B'(z"?dz"dz' . 



(3.51) 



^0 
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where we used ^(0) = to remove the integration constant. As -B(l) diverges at 
the extremal hmit, B{z) would diverge logarithmically as i? ~ — ln(l — z). This 
implies that the r.h.s. of Eq. (3.51) also negatively diverges at least logarithmically, 
and then A{1) -> -oo, or e^^^^ 0. 

4. Einstein- Maxwell-dilaton-axion SL(2,R) model 

In this section, we consider Einstein-Maxwell theory coupled to a dilaton-axion with 
the action 

S = j d^x^ i^^lj' 2( V0)' - ^e^'^(Va)2 - e'^'^F^ - aFF^ . (4.1) 
It is well known that this action is invariant under global SL{2,R) transformations 



A -> , dd-bc=l, (4.2) 

where A = a + ie~^'^ = Ai + 1X2 and 

F^, ^ (cAi + d) F^, - cX^F^, . (4.3) 

Therefore, once we obtain a charged black brane solution of the theory, we also obtain 
both electrically and magnetically charged black brane solutions by the SL{2, R) 
transformations. This model is studied in holographic setting for the Quantum Hall 
effects in [27, 41]. At the quantum level, this SL(2,R) symmetry will be enhanced 
to SL{2,Z). For the connection of SL{2,Z) to Quantum Hall Effects, see [42, 43, 
44, 45, 46] and review [47]. 

In the purely electrically charged case, the last term in Eq. (4.1) vanishes and 
then, the equations of motion are 

R^u = -^9>.u + 2V^(t>VA + ^e^-^V^aV.a + 2e-^^ F,,^F,'' - ^g^^e'^'^F^ , (4.4) 

□0 _ le^<^(Va)2 + le-2<^F2 = , (4.5) 

Da + 4V;,0V^a = , (4.6) 
V^V.ie-^'^Fn = d.iV^e^'^Fn = . (4.7) 

We take an ansatz 

= 0(2), a = ax, A^dx^ = A^dv, c{z) = , (4.8) 

under the metric (3.5). Eq. (4.6) is automatically satisfied and the solution of 
Eq. (4.7) is written by 

F,, = a; = -C2Le- V^ (4.9) 
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where C2 is a constant corresponding to the charge density of the gauge field. The 
energy- momentum tensor is 

-^5'', (^2(V0)' + ^e^<^(Va)2 + e'^^F^^ . (4.10) 

Then, ansatz (4.8) satisfies T^^ T^x = O5 therefore we expect to have a 

B{z) 7^ and c{z) = solution. 

Substituting Eq. (4.9) into Eqs. (4.4) and (4.5), we obtain four-coupled differen- 
tial equations, 

2A" + A'^ + B'^ + 4(j)'^ = , (4.11) 

zgB" + {zg' + g{zA' - 2)}B' + ^a^e^'^-^^'' z = , (4.12) 

2^(2 - zA')g' - {12 - 8zA' + z\A'^ - B'^ - A(j)'^)}g 

+12 - «2;j2g40-A-B _ 4c2/e2<^e-2^ = , (4.13) 

e^g<P" + z' (^^^ ' 0' - ^a'e^'^e-^ - clz'e-^e'^ = , (4.14) 

where the third equation corresponds to the constraint equation. The equations of 
motion (4.11) - (4.14) are invariant under the transformation 

^ - 00 , a ^ e^-^^a , C2 ^ e'^Ocs , (4.15) 

for an arbitrary value of 0o. This is nothing but the SL{2, R) symmetry given by 
(4.2) and (4.3) with SL{2, R) parameter a = = e'^°, b = c = 0. Using the freedom 
of this transformation (4.15) and the rescaling of the coordinate x and y, we can set 

A(l)=i?(l) = 0(l) = (7(l) = O. (4.16) 

By Eqs. (4.11) - (4.14), the regularity of the black brane solution requires 

Am = 2^'l^^, B'(i) = |^. m^-'-^, (4.17) 

where g'{l) = —n- Thus, the regular solutions of Eqs. (4.11) - (4.14) are uniquely 
determined by three parameters, a, C2, and k. We consider the boundary condition 
(3.46) adopted in section 3. For simplicity, we also require that = at the infinity. 
Since the asymptotic behavior of B is given by Eq. (3.49), we can set A{0) = B{0) = 
-B'(O) = by the coordinate transformation (3.48). So, we can numerically find the 
value of C2 satisfying A'{0) = at the infinity for a given k and a. Once the solutions 
are obtained, we can set 0(0) = by the transformation (4.15). 

Figs. 7-12 show the anisotropic black brane solutions for the rescaled coordi- 
nates. Figs. 7-9 plot the horizon area e"^^^-*, and the values of B, at the horizon 
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Figure 7: (color online) K-e"^^^) relation. The area of the horizon per unit Ax = Ay = 1 
are shown as a function k for a = The area rapidly decreases to zero as k — )• 0. 
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Figure 8: (color online) k-B{\) relation for a = \f2. B{1) goes to zero as k decreases to 
zero. 
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Figure 9: (color online) K-(f){l) relation for a = \/2. 
to zero. 



(1) negatively diverges as k decreases 



as a function of k (the minimum of k is 0.01) for a = \/2. Due to the existence of 
the dilaton field (p, the horizon area rapidly decreases to zero at the extremal limit, 
K — )■ (Fig. 7), and the anisotropy also decreases (Fig. 8.), contrary to the previous 



- 19 - 



D 




Figure 10: (color online) D = z^g^xj^^^ (blue), z^gyyjl? (green), and (red) are shown 
as a function of z for a = \/2) k = 1. 
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Figure 11: (color online) D = z'^g^x/L'^ (blue), z'^gyy/L'^ (green), and (j) (red) are shown 
as a function of z for a = \/2, k = 0.1. 
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Figure 12: (color online) D = z'^gxx/L'^ (blue), z'^gyy/L?' (green), and (j) (red) are shown 
as a function of z for a = \/2, k = 0.01. 



real scalar model. This is because the negative divergence of the dilaton field strongly 
suppresses the kinetic term e^''^(Va)^ in Eq. (4.1) which generates the anisotropy. 

Figs. 10 - 12 show the metric as a function of z for k = 1, 0.1, 0.01 and a = ^J2. 
We can see that the difference between g^x and gyy goes to zero as k decreases to 
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zero. This means that the geometry near the horizon becomes isotropic when the 
black hole approaches the extremal solution. 

Given the solution shown in Figs. 7 - 12, then by SL{2, R) symmetry, we can ob- 
tain more generic solutions where all the fields, i.e., dilaton 0, axion a, field strength 
Fzv and F^y have more generic position x-dependence. However, since SL{2, R) sym- 
metry does not transform the metric, the metric is still anisotropic but homogenous. 
This means that the position dependence of 0, a, F^v, F^y does not yield the position 
dependence for the energy- momentum tensor. 

5. The third law of anisotropic black branes 

In this section, we consider the third law of thermodynamics for the general ain- 
sotropic black brane solutions. The third law states that the area of the black brane 
horizon goes to zero as the temperature decreases to zero. As seen in Eq. (3.51), the 
area must go to zero when B{z) diverges. Indeed, we can give a theorem supporting 
the third law of thermodynamics. 

We consider the following static anisotropic black brane solutions with the metric 

ds^ = ^ (^~giz)dv^ - 2dvdz + e^(^)+^(^)rfx2 + e^^^^-^^^^rfy^) . (5.1) 
As the asymptotic boundary condition, we shall impose that 

lim A{z) = lim B{z) = lim A'{z) = lim B'{z) = 0. (5.2) 

2—^0 2— 5>0 2— >0 Z~^0 

In the neighbourhood of the horizon, we assume that there is a small positive 
value e such that inl — e<2;<l,5f can be expanded as 

g{z) = k{1 -z)+ g2{K){l - zf + 0(6=^). (5.3) 

Note that on top of the horizon where 2^-1, the first term always dominates over 
the second term. However, slightly away from the horizon, like z = 1 — e, with e > 0, 
then, near the extremal limit where k — )■ 0, the second term dominates over the first 
term. Given this, then we can show the following theorem: 

Theorem. 

For the anisotropic black brane solutions with metric (5.1) satisfying the Einstein 
equations, if the following two conditions, 

1. The null convergence condition is satisfied, i.e., R^uV^V^ ^ for any null 
vector V^, 

2. Defining 

Ar(.) = f {Ty^{z) - T\{z)) = e-^+%y - e-^-^T.. (5.4) 
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as the difference for the re-normahzed^ energy-momentum-tensor, it is ex- 
panded as 

AT{z) = a{K) + -f{K){l - z) + 0{e^) (5.5) 
inl-e<2;<l and hm«^+o a^n) = ao (7^ 0), \im^^+o-f{K) = 70, 

are satisfied, the horizon area must go to zero as k decreases to zero. 

The second condition imphes that anisotropic matters, which induces nonzero 
B{z), exist near the horizon. We shall prove the theorem by using the following 
Lemma. 

Lemma. 

A is non-positive and A and u = e"^ satisfy the following inequalities, 

u' 

< u{z) < 1, e = - = A' <0. (5.6) 

u 

proof) 

From the Einstein equations, we obtain 

e' = -^e'-T,,-^B". (5.7) 

From the condition 1, we have T^z = Rzz > 0. Thus, A' must be a non-increasing 
function. By the asymptotic boundary condition (5.2), i.e., A{Q) = A'{0) = 0, A' < 
and A < 0. In addition, by definition of u, we have < u{z) < 1. □. 

proof of Theorem). 

Let us suppose that the area at the horizon does not go to zero at the extremal hmit, 

lim A{1; k) = Aq. (5.8) 

Integrating (2.24) by z, we obtain 

where the constant of integration is determined by imposing regularity at the horizon. 



/ is due to the overall metric factor 
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When K is small enough, by the condition 2, the sign of ^T{z) does not change 
for 1 — e < 2; < 1. So, \B'iz)\ can be evaluated for 1 — e < z < 1 as 



\B'{z) 



z e 



^ ATjz') 
\AT{z') 



y/2 



e^^'^dz' 



> 



> 



z e' 



-Aiz) 



' \AT{z') 
^12 



,2gA(l)-A(l-e) rl 



yl2 



-dz' 



Uz'. 



(5.10) 



Here, we used Lemma for the first and second inequalities. Substituting Eqs. (5.3) 
and (5.5) into (5.10), Eq. (5.10) is reduced to 

u{l){\a{K)\ + 0{e)) 



\B'(z)\ > 



(5.11) 



u{l - e)iK + g2iK){l - z)) 

for 1— e < 2 < 1. Note that at the extremal limit k — 0, we can make \a{K,)\/{K, + g2{n){l — z)) 
as big as 1/k as we go near the horizon z — )■ 1. This indicates that B'{z) diverges as 
1/(1 — z) and B{z) diverges logarithmically near the horizon in the extremal limit. 
One exception for this is, u{l) — >■ 0. However, we have assumed this is not that case 
by (5.8). 

Since Tzz = Rzz > by the condition 1, Eq. (5.7) means 

1 



0' < 
- 2 



(5.12) 



Therefore, B' diverges as 1/(1 — z) induces that 6 also diverges negatively as — 1/(1 — 
z). This indicates that A diverges logarithmically to negative infinity in the extremal 
limit, therefore — )■ 0, but this contradicts with (5.8). 

To say this in more rigidly, using Eq. (5.11), we can integrate Eq. (5.12) from 
u=l— etou = l and thus we obtain 



9{z)-9{l-e) < 



2^2(1 — e)ng2{K) 



1 



1 



32 (ft) I 



K J 



(5.13) 



Since 6'(1 — e) < by Lemma, it immediately means that 



iz) < 



u\l){a\K) + 0{e)) 



1 + ^(1 



2^2(1 — e)K,g2{K,) 
Integrating Eq. (5.14) again from u = 1 — e to u 



Z) 1 + 

1, we finally obtain 



(5.14) 



n{l) 

m2(1) "^^(i-e) 



In 



< 



a^{K) + 0{e) 
2^2(«:) 



.92{l^) 



In 1 + 



92[,i^)e 



K 



K + eg2{K) 



(5.15) 
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This indicates that there is a small positive value such that 



u2(l-e), u(l) 

^— ^In— ^<-D, 0<K<Ko 5.16 

for an arbitrary large positive value D. By assumption (5.8), ^ for small 
positive value n = ki < hq. This would be u{l — e; ni) > 1. This contradicts with 
the fact that u(l — e) < 1 in Lemma. □. 

Therefore, at the extremal limit, (5.8) gives the contradiction and we conclude 
that area should approach zero at that limit. This is consistent with the numerical 
solutions we obtained in section 3 and 4, Fig. 1 and 7. 



6. Discussion 

In this paper, we have studied several four dimensional gravitational theories and we 
have obtained homogeneous but anisotropic black brane solutions given by the metric 
ansatz (5.1), which asymptotic to AdS^ space-time. The gravitational systems we 
consider are 1) Einstein-Maxwell dilaton theory, and 2) Einstein-Maxwell-dilaton- 
axion theory with SL{2, R) symmetry. The anisotropy is induced in our setting by 
the scalar field profile, either by dilaton or axion a. The scalar fields have manifest 
origin of anisotropy and induce the anisotropy for the metric while keeping the metric 
homogeneous at fixed radial slice. However, the scalar fields themselves have manifest 
inhomogeneous profile. In the case 1), we obtain solutions both analytically and 
numerically. Analyitic solutions are obtained by the perturbation from the isotropic 
RN black brane solution, and in the case 2), we obtain solutions numerically. Our 
solutions are smooth everywhere but seem to approach singular behavior at the 
horizon in the extremal limit, as seen from Fig. 6, 12. 

We also showed how the third law of thermodynamics holds in our set-up, which 
implies that as temperature goes to zero, the area also becomes zero. Strictly speak- 
ing, at the limit temperature goes to zero, the solution seems to approach singular 
behavior. As far as we have analyzed, we could not obtain the regular black brane 
solution analytically in the extremal limit. It is the numerical analysis which suggests 
that, the solution approaches singular in the extremal limit. 

Furthermore, even if we do not take the zero temperature limit, as we lower the 
temperature, the curvature of the black brane solutions become bigger and bigger 
near the horizon, and there is a critical temperature where two- derivative Einstein- 
Hilbert action breaks down. Then, we have to worry about higher derivative correc- 
tions coming from stringy corrections or quantum corrections. Therefore, precisely 
speaking, the extremal limit may not be meaningful limit in our two derivative Ein- 
stein Hilbert action. This suggests that we may always need non-zero but small 
temperature as an IR regulator for the system. 



-24- 



Because we have used the Einstein equations in our proof of third law in section 
5, our argument at that section is vahd for two-derivative Einstein-Hilbert action. In 
holography, this two-derivative assumption in bulk is a big assumption which holds 
only the strong coupling limit in boundary field theory. If the coupling constant in 
field theory is not that big, strong coupling expansion in field theory side always 
induces stringy corrections in the gravity side. Therefore we should rather regard 
our argument of third law of thermodynamics as an "indication" that as we lower 
the temperature more, the horizon area becomes smaller and smaller. It is curious 
to develop our argument furthermore without relying on the two-derivative action 
assumption. 

There are several related questions remaining for such anisotropic solutions. One 
of those is checking the stability. It can happen that such anisotropic solutions be 
unstable by the small perturbations. If this is the case, then such solutions give at 
most meta-stable phases in the dual field theory. We left this for future study. 

In our paper we have restricted our analysis only to the metric of the form 
(5.1), namely we have set c{z) = in more generic metric (2.1). This is because we 
are considering the scalar fields which give T^^ = 0. Obviously with more generic 
matter contents, we can also obtain anisotropic solutions which have nonzero c{z). 
Generic solutions, for examples, as the ones analyzed in the Bianchi type [33] does 
not always satisfy c{z) = 0. Therefore, it is interesting to search more generic setting 
with c{z) ^ 0. 

As is discussed in section 2, this implies that we need the matter contents which 
satisfy T^^ 7^ in addition to T^y 7^ T^^. For example, if we introduce a scalar field 
profile such that 



we can obtain both T^^ 7^ and T^y 7^ T^^. However, generically such scalar fields 
induce inhomogeneous metric, this is because energy-momentum tensor T'^^ becomes 
manifestly position-dependent. By requiring that an energy-momentum tensor is 
position- independent, and that it satisfies T^^ 7^ and T^y 7^ T^^, we end up with 
a scalar field profile such that it can depend on coordinate x and y only linearly. 
But this is, after coordinate transformation, equivalent to our dilaton profile (3.9). 
Therefore if we introduce anisotropy by the scalar fields, our setting is quite generic 
one in that sense. 

Of course, we can also induce anisotropy from other degrees of freedom, such 
as the gauge potential [18], [33], and generical p-foim potential [39]. Generically, 
the gauge potentials or p-form potentials which are invariant under Killing vectors 
respect the isometry of the metric and therefore do not induce inhomogeneity. Taking 
the metric and gauge potential made by the one forms of the Bianchi classification, 
we can obtain homogeneous but anisotropy solutions as [33]. However we can also 



= f{x,y) , 
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take the gauge potential, say Ay to be proportional to x 

Ay = ax , (6.2) 

as we took for the scalar fields in this paper. In such case, the field strength is 
position x- independent, so is the energy-momentum tensor. Therefore, we can obtain 
homogeneous but anisotropic metric. However, even though we have homogeneous 
metric, the gauge potential is manifestly x-dependent and we have a source term, 
i.e., constant flux on the boundary theory. We can also construct homogeneous 
but anisotropic solutions with the scalar flelds and vector flelds where scalar flelds 
and gauge potentials are not simply proportional to x. By acting the SL{2, R) 
symmetry (4.2) and (4.3) on the solutions we obtained in section 4, we can have 
generic solutions where scalar flelds and gauge potentials have non-trivial position 
x-dependence. Obviously, it is desirable to investigate more generic homogeneous 
and anisotropic solutions where anisotropy is induced not only by the scalars but 
also by the generic gauge or p-form potentials in more detail. 

It would be also interesting if we could construct the situation where anisotropy 
is induced spontaneously as [39] , instead of the situations where anisotropy is induced 
by the source terms (non-normalizable modes) such as (6.1) and (6.2). 

The metric we consider is such that, the fleld theory spaces, spanned by x and 
y, becomes rectangular, namely the ratio between x and y directions changes due to 
the nonzero B{z). This is the situation where one direction is stretched compared 
with the other. For example, if we consider the fermion Green function, then the 
Fermi-surface should change into elliptical shape instead of circle in momentum {kx, 
ky) space. Since our solutions constructed in section 3 and 4 have non-trivial gauge 
potentials A^ which is dual to the chemical potential in boundary theory, it is quite 
interesting to calculate the fermion Green function in such anisotropic set-up as 
[48, 49, 50, 30]. It is also interesting to study the transport coefficients in such 
anisotropic settings. We left these studies as future projects. 
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A. Analysis of uniqueness inside the event horizon 

In section 2.2, we have shown that i? = c = 0in2;<lin the pure gravity system. 
This can be extended just inside the horizon as follows. Let us suppose that both B 
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and c are not zero but small just inside the horizon, < z — 1 <^ 1. Then, Eqs. (2.6) 
and (2.7) are reduced to the following linear second order differential equation, 

gK" + g'K' = 0, (A.l) 

where K = B, c and we used the fact that g, c are negligible with respect to g', in 
< z — 1 <^ 1. Since g ~ k(1 — z) near the horizon, Eq. (A.f ) has a regular singular 
point at 2; = 1. Thus, substituting an ansatz K = {1 — zf into Eq. (A.f), we obtain 
7^ = 0. This means that the two independent leading solutions are 

~ f , K2^\^\\-z\. (A.2) 

Using two arbitrary constants Ci, C2, the general solution is expressed as ~ 
CxKx + C2K2 in < z - f < f. As = 5, c = at z = f , we must take 
C\ = C2 = 0. This indicates that anisotropy does not still appear even just inside 
the horizon, i. e. , 5, c ~ when 2; — f -C f . 



B. Scalar curvature invariants for the generalized Gaussian 
null coordinates 

Given the generalized gaussian null coordinate (2.f ), the scalar curvature is calculated 
as 



2L2 (e2^W - c{zYf 

where 

Kr{z) = Az^e^^'^'h{zfg{z)A"{z) - Az^e''^^''>g{z)A"{z) - 8z^e^^^'^c{z)g{z)A'{z)c'{z) 
+Az^e^^^'^cizfA'{z)g'{z) + 7z^e'^^'h{zfg{z)A\zf - 12ze^^^'^cizfg{z)A'{z) 
-Az^e^^^'^A'{z)g'{z) - 3z^e''^'^''>g{z)A'{zf + 12ze^^^'^g{z)A'{z) 
+z^e^^^'^c{zyg{z)B'{zf - z^e^^^'^g{z)B'{zy + Az^e^^^'^c{z)g{z)c"iz) 
+Az^e^^^''^c{z)c'{z)g'{z)+3z^e^'^^''^g{z)c'{zf - 12ze^^^'^c{z)g{z)c'{z) 
^^z^e^^^^^c{zfg"{z) - 2Aze^^^''>c{zfg'{z)+A8e^^^''>c{zfg{z) - 2z^e^^^'^g"{z) 
+12ze^^^'^g'{z) - 2Ae''^^'^ g{z) - Az^c{zfg{z)d\z) - Az^c{zfd {z)g\z) 
+z^c{zfg{z)c'{zf + 12zc{zfg{z)c'{z) - 2z^c{zYg"{z) 

+ l2zc{zYg\z) - 2Ac{zYg{z) . (B.2) 

Note that the curvature is non-divergent if Kr{z) is non-divergent and 
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The former condition is satisfied if A, A', A", c, c', c", g, g', g", B' are regular 
functions. The latter condition implies that the determinant of the matrix gij = 
X, y) on the v = const, two dimensional spacelike surface is positive (the fixed v slice 
has positive area). 

Similarly, Ricci tensor square R^'^R^y and the Kretschmann scalar curvature 
invariant R^"'''^R^yap are calculated as 

and Kfiicci{z) and Kjnemanniz) are similar non linear combination of A, A', A", c, d 
c", g^ g\ g'\ B', B", and non-divergent if these are so. Therefore, with (B.3) if the 
functions A{z), B{z), c{z) g{z) are all smooth up to their second derivatives, all the 
scalar curvature invariants are finite. 
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